function ny = func_ny_from_nxy(nz, nx, CONSTS)

    eps_a = CONSTS.eps_a;

    ny = sqrt_from_negative_num(eps_a^2-nx.^2-nz.^2);

end